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THE  MISSION  OF  AGARD 


The  mission  of  AGARD  is  to  bring  together  the  leading  personalities  of  the  NATO  nations  in  the  fields  of 
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-  Continuously  stimulating  advances  in  the  aerospace  sciences  relevant  to  strengthening  the  common  defence 
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-  Improving  the  co-operation  among  member  nations  in  aerospace  research  and  development; 
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field  of  aerospace  research  and  development; 

Rendering  scientific  and  technical  assistance,  as  requested,  to  other  NATO  bodies  and  to  member  nations 
in  connection  with  research  and  development  problems  in  the  aerospace  field, 

-  Providing  assistance  to  member  nations  for  the  purpose  of  increasing  theii  „_ientific  and  technical  potential; 

-  Recommending  effective  ways  for  the  member  nations  to  use  their  research  and  development  capabilities 
for  the  common  benefit  of  the  NATO  community. 

The  highest  authority  within  AGARD  is  the  National  Delegates  Board  consisting  of  officially  appointed  senior 
representatives  from  each  member  nation.  The  mission  of  AGARD  is  carried  out  hrough  the  Panels  which  are 
composed  of  experts  appointed  by  the  National  Delegates,  the  Consultant  and  Exc  lange  Program  and  the  Aerospace 
Applications  Studies  Program  The  results  of  AGARD  work  are  .eported  to  the  member  nations  and  the  NATO 
Authorities  through  the  AGARD  series  of  publications  of  which  this  is  one. 

Participation  in  AGARD  activities  is  by  invitation  only  and  is  normally  limited  to  citirer.s  of  the  NATO  nations. 


The  material  in  ;his  publication  has  been  reproduced 
directly  from  ccpy  supplied  by  AGARD  or  the  author. 


Published  December  1973 


620.178.3:539.431:624.07 


$> 

Printed  by  Technical  Editing  and  Reproduction  Ltd 
Harford  House.  7  ■  9  Charlotte  St  London.  W IP  I  HD 


This  volume,  the  third  part  of  a  series  giving  data  for  design  against  acoustic  fatigue, 
has  been  prepared  in  order  to  draw  together  the  results  of  research  in  acoustic  fatigue  and 
to  present  them  in  a  form  directly  usea'ole  in  aerospace  design.  Future  work  in  this  series 
will  deal  with  endurance  of  titanium  alloy  structures  under  simulated  acoustic  loading,  near 
field  compressor  noise  estimation,  stress  response  of  box  and  control  surface  structures, 
structural  damping  and  stress  response  of  skin-stringer  panels  with  stringers  of  relatively 
low  flexural  stiffness. 

T  he  AGARD  Structures  and  Materials  Panel  has  for  many  years  been  active  in 
encouraging  and  coordinating  the  work  that  has  been  necessary  to  make  this  collection 
of  design  data  possible  and  after  agreeing  on  procedures  for  the.  acquisition,  analysis  and 
interpretation  of  the  requisite  data,  work  on  this  series  of  design  data  sheets  was  initiated 
in  1970. 

The  overall  management  of  the  project  has  been  conducted  by  the  Working  Group 
on  Acoustic  Fatigue  of  the  AGARD  Structures  and  Materials  Panel,  and  the  project  has 
teen  financed  through  a  collective  fund  established  by  the  Nations  collaborating  m  the 
project,  namely  Canada,  France,  Germany,  Italy,  UK  and  US.  National  Coordinators 
appo...ted  by  each  country  have  provided  the  basic  data,  1L  ised  with  the  sources  of  the 
data,  and  provided  constructive  comment  on  draft  data  sheets.  These  Coordinators  are 
Dr  G.M.Lindberg  (Canada),  Mr  R.Loubet  (France),  Mr  G.Bayerd&rfer  (Germany), 
Gen.A.Griselli  (Italy),  Mr  N.A.Townoend  (UK),  Mr  A.W.Kolb  (US)  and  Mr  F.F.Rudder 
(US).  Staff  of  the  Engineering  Sciences  Data  Unit  Ltd,  London,  have  analysed  the  basic 
data  and  prepared  and  edited  the  resultant  data  sheets  with  invaluable  guidance  and 
ad  'ice  from  the  National  Coordinators  and  from  the  Acoustic  Fatigue  Panel  of  the 
Ro/al  Aeronautical  Society  whicn  has  the  following  constitution:  Professor  B.L.Clarkson 
(Chairman),  Mr  D.C.G. Eaton,  Mr  J.A.Hay,  Mr  W.T.Ktrkby,  Mr  M.J.T.Smith  and 
Mr  N.A.Townsend.  The  members  of  staff  of  the  Engineering  Sciences  Data  Unit 
concerned  with  the  preparation  of  the  data  sheets  in  this  volume  are:  Mr  A.G.R .Thomson 
(Executive,  Environmental  Projects),  Dr  G.Sen  Gupta  and  Mr  R.F. Lambert  (Environmental 
Projects  Group). 

Data  sheets  based  on  this  AGARDograph  will  subsequently  be  issued  in  the  Fatigue 
Series  of  Engineering  Sciences  Data  issued  by  ESDU  Ltd,  where  additions  and  amend¬ 
ments  will  be  made  to  maintain  their  ament  applicability 


A.ri.Haii 

Chairman, 

Working  Group  on  Acoustic  Fatigue 
Structures  and  Materials  Panel 
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Section  1 

ENDURANCE  OF  TITANIUM  AND  TITANIUM  ALLOY  STRUCTURAL  EJJSMENTS 
SUBJECTED  TO  SIMULATED  ACOUSTIC  LOADING 

1 . 1  Notati an 

S  root  mean  square  value  of  street  at  a  reference 

r**  position 

N  equivalent  endurance,  taken  as  half  the  number 

r  of  Kero  crossings  of  the  stress-time  function 

before  failure 

1.2  Notes 

This  Section  gives  the  results  of  fatigue  tests  on  titanium  alloy  specimens, 
typical  of  aircraft  structural  elements,  excited  by  narrow-band  loading  of  random 
amplitude  with  zero  moan  load  to  t Insulate  stress  response  to  acoustic  loading. 

The  plain  specimens,  illustrated  in  Figure  1.1a,  are  representative  of  sheet 
material  away  from  the  effects  of  attachments.  The  built-up  specimens,  illustrated 
in  Figures  1.1b  to  l.le  represent  conditions  in  sheet  material  attached  without 
Joint;. ng  compound  to  stiffeners,  where  adjacent  plates  are  vibrating  in  phase. 

In  Figures  1.2  to  1.5,  S  ie  plotted  against  N  for  teat  results  grouped 

rms  r 

according  to  the  type  of  structural  element,  material  and  method  of  attachment  as 
shown  in  Table  1,1.  A  etraight  line  'itted  by  the  least  squares  method  and  bands 
enclosing  approximately  95  per  cent  oi  the  data,  (two  standard  deviations  either 
side  of  the  mea/i'  have  been  given  for  guidance;  they  should  not  be  used  directly 
<"or  estimating  life  because  due  account  oust  be  taken  of  scatter.  Figure  1,6  shows 
tne  relative  endurances  for  the  different  groups.  Extrapolation  of  the  mean  lines 
beyond  the  limits  Indicated  on  Figure  1.6,  where  few  data  arc  available,  is  expected 
to  lead  to  an  underestimate  of  endurance.  Data  from  similar  tests  under  constant 
amplitude  loading  indicate  that,  for  titanium  alloy3,  the  S-N  curve  tends  to  level 
7 

off  above  10  cycles. 

The  reference  position  for  sriro  taken  as  a  point  on  the  failure  line  away 

from  the  effocts  of  stress  concentration  or  fxet*ing  at  fastenings;  the  test  points 
plotted  represent  either  mean  stresses  over  the  areas  r  . verod  by  monitoring  strain 
gauges  affixed  across  the  failure  line  away  from  fast-  -.ings,  or  stresses  at  that 
position  estimated  after  calibration  from  .strain  monl.ors  located  'sewhere.  When 
uolng  this  Section  to  estimate  a  life  using  a  root  a  .'an  square  stress  level  obtained 
from  Reference  1.4.5,  It  is  recommended  that  the  calculated  ’-alue  of  stress  at  the 
rivet  line  should  be  assumed  to  bo  the  same  as  the  stress  at  the  failure  line.  In 
practice  tho  failure  position  is  sufficiently  cluse  to  the  rivet  lino  'or  this 
approximation  to  be  within  the  range  of  accuracy  of  the  simple  theory  used  for  stress 
prediction  in  Reference  1.4.5. 

Details  of  the  test  methods  and  specimens  are  given  in  appendices  as  listed 
below. 


N/m2  lbf/in2 
cycles  cycles 


Appendix  1A.1  -  Tests  on  plain  specimens 
Appendix  1A.2  -  Tests  on  fa3tened-skin  specimens 
Appendix  1A.3  -  Di  scnptlon  of  materials 

Typical  failure  locations  are  shown  on  Figure  1.1. 

The  loading  simulated  fc  *■  the  fastened-skl n  specimens  is  only  that  for  modes 
where  adjacent  panels  are  in  phase.  Care  should  be  taken  in  using  the  data  for  modes 
involving  twisting  of  the  stiffeners. 

1.3  Interpretation  Of  The  Results 

1.3.1  Effect  of  method  of  attachment 


Skin  specimens  attached  to  stiffeners  by  mushroom-headed  rivets  showed 
significantly  grsater  snduranoe  and  less  scatter  than  resistance  spot-welded 
specimens.  The  resistance  spot-welded  specimens  in  turn  showed  slightly 
greater  endurance  and  considerably  less  sertter  than  gas  tungsten-arc  spot-woldod* 
specimens.  The  latter  pro res  therefore  requires  furthor  development  before  it 
can  be  used  with  confidence  for  titanium  structures  in  a  high  noise  environment. 


*  The  gas  tungsten-arc  spot  welling  process  is  described  in  Reference  1.4.4. 
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1.3.2  Effect  of  material  s  pacification  and  thickness 

The  available  data  for  different  materials  are  not  direotly  comparable 
because  of  difference#  in  specimen  configuration  and  thickness.  The  large 
differences  betveon  the  result#  for  Ti-2Cu  (Aged)  (Figures  1.3  and  1,4)  and 
C.P.  Ti  (Figure  1.5)  are  thought  to  be  due  <  ->  the  effeots  of  skin  thickness 
on  spot-weld  performance  rather  than  to  the  difference  in  materials, 

1.3.3  Effect  of  heat  treatment  and  grain  direction 

The  tests  on  the  Ti-2Cu  material  show  that  for  this  material  a  substan¬ 
tial  increase  in  endurance  is  obtained  by  spot  welding  before,  rather  than 
after,  ageing.  For  those  specimens  slightly  greater  endurance  was  obtained  for 
spot  welds  stressed  in  a  transverse  rather  than  a  longitudinal  grain  direotion 
but  the  difference  might  not  be  significant. 

1.3.4  Effeot  of  test  frequency 

The  mean  responee  frequency  for  eaoh  type  of  test  specimen  is  indicated 
in  Figures  1.2  to  1.4.  Tha  range  of  frequencies  is  too  small  for  any  effects 
to  be  identified. 

1.3.5  Effect  of  elevated  temperature 

A  limited  number  of  tesis  on  Ti-2Cu  (Aged)  specimens  was  carried  out  at 

250°C.  These  indicate  that  temperatures  up  to  250°C  have  little  effect  on  the 
endurance  of  this  material. 

1.3.6  Data  for  plain  speelnsns 

Tests  on  the  plain  specimens,  under  constant  amplitude  loading,  showed 
a  very  marked  discontinuity  in  the  S-N  curve  at  alternating  stresses  of  about 

2  2 

450  MN/m  (65  000  lbf/in  ).  The  increased  scatter  of  the  test  results  at 

2  2 

r.m.a.  stresses  of  about  250  MN/m  (35  000  lbf/in  )  in  Figure  1.2  is  believed 
to  be  associated  with  this  discontinuity. 


1.4  Derivation  and  References 
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TABLE  1.1 

SUMMARY  OF  DATA  PRESENTED 


Figure 

No. 

Type  of 
E'ement 

— 

Material 

Number  of 
Attachment 
Rowa 

Attachment 

Type 

Thickne»3 

Response 

Frequency 

<H«) 

mm 

in 

1.2 

Plain 

T1-6A1  4V 

1.52 

(0.0 6) 

300 

(Ann) 

.  _ 

1.3 

Riveted 

1 

c's’k  Monol 

(0.71) 

0.028 

260 

skin 

Ti-2Cu( Aged) 

rivet 

2 

nt/h  Monel 

(0.71) 

0.028 

160 

rivet 

1.4 

Welded 

1 

resistance 

(0.71) 

0.028 

250 

skin 

spot  weld 

1.5 

Welded 

C.P.Ti 

2 

resistance 

(0.41) 

0.016 

190 

skin 

spot  weld 

2 

gas 

(0.41) 

0.016 

180 

1 

spot  weld 

(0.41) 

0.016 

175 

1.6 

ASSEMBLAGE  OF  DATA  PLOTTED  ON  FIGURES  J 

.2  TO 

1.5 
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1A, 1  T»gt»  On  Plain  Sp#oi»«ng 
Method  of  leatiaj: 


APPENDIX  1A 


Plain  cantilever  specimens  of  the  type  shown  In  Figure  1.1a,  rigidly  clasped  at 
one  end,  were  excited  electromagnstloally  through  small  ateel  platee  attached  to  the 
other  ends  auch  that  they  vibrated  at  their  resonant  frequency  of  approximately 
300  Hz  with  a  Rayleigh  distribution  of  peak  stresses  truncated  at  ±4  tines  the  root 
mean  square  stress.  The  stress  level  in  the  specimens  was  monitored  by  a  capacitive 
transducer,  located  about  half-way  along  the  specimen,  which  had  been  calibrated 
against  strain  gauges  affixed  at  the  failure  location. 


The  random  loading  was  interrupted  at  intervals  to  test  for  failure,  using 
constant  amplitude  excitation,  the  criterion  being  the  reduotion  of  the  resonant 
frequency  to  98£  of  its  original  value. 


1A.2  Teats  On  Fastened -Skin  Specimens 

Method  of  testing 

Skin-stiffener  specimens  of  the  types  shown  in  Figure  2  were  clamped  by  the 
base  to  a  moving  coil  vibrator  fed  from  a  white  noise  generator.  The  signal  from 
the  noise  generator  was  filtered  to  a  $  octave  bandwidth  centred  at  tho  fundamental 
resonant  frequency  of  the  specimen,  so  giving  rise  to  vibration  of  the  speoimen  with 
random  amplitude  at  its  resonant  frequency.  The  truncation  stress  level  was  between 
±3.0  and  ±3.3  times  the  root  mean  square  level.  The  stress  level  in  the  specimens. 
was  monitored  during  the  tests  by  strain  gauges  located  either  on  the  expected 
failure  line  or  on  the  centre-line  of  the  fasteners.  In  the  latter  case,  failure-line 
stresses  were  estiaated  using  a  factor  determined  from  measurement  of  dynamic  strain 
distribution  across  the  Joint. 


Failure  was  detected  by  the  increase  in  damping  and  a  marked  change  in  response 
characteristics . 


1A.3  Description  Of  Materials 
General  descriptions 

The  designations  used  in  this  Section  are  descriptive  of  the  chemical  composition 
and  condition  of  heat  treatment  of  the  material. 


C.P.Ti 

Commercially  pure  titanium.  The  material  used  In  the  tests  was  the  high  strength 

grade  (tensile  strength  620  to  770  MN/n»2  (89  600  to  112  000  lbf/in2))  with  commercial 
designation  IMi  l60.  The  UK  specification  for  this  material  was  D.T.D.  5063,  but  this 
specifioation  was  supersedad  in  1968  by  D.S.  TA6  whioh  covers  a  slightly  lower  range  of 
tensile  strength.  This  material  is  similar  to  ths  US  commercial ly  pure  titanium, 

T1-75A  (Specifioation  MIL-T-9046E  Type  1  Comp.B  or  AMS  4901). 

Tl-2Cu  (Aged) 

This  is  a  weldable  binary  titanium  alloy,  susceptible  to  ags  hardening, 
containing  2  to  3  per  cent  copper.  The  specification  for  the  material  used  in  t».e 
tests  was  BAC  M.47  and  the  commercial  designation  was  IMI  230  (Aged).  The  design 

tensile  strength  for  this  material  is  695  MN/m2  ( 100  800  lbf/in2).  The  UK  specifioation 
for  this  type  of  material  is  D.T.D,  5233.  There  is  no  US  Military  or  AMS  specification 
covering  a  similar  material. 

T1-6A4  4V  (Ann) 

The  material  tested  had  the  nominal  ohomloal  composition:  6  per  cent  aluminium, 

4  per  cent  vanadium,  0.08  per  cent  carbon,  low  nitrogen  and  hydrogen,  remainder 

2 

titanium.  The  meohanioal  properties  were.  tensile  strength,  960  MN/m 

(133  000  lbf/in2)!  yield  strength,  890  KN/m2  (129  000  lbf/in2) |  elongation, 

10  per  cent. 

The  material  is  similar  to  US  apecif ioation  MII-T-9046E,  Type  III,  Comp.C 
(Annealed)  or  AMS  4911  (Annealed)  or  British  Standard  TA  10. 
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Section  2 

DAMPING  IN  ACOUSTICALLY  EXCITED  STRUCTURES 


2 . 1  Note* 


Duping  is  an  important  parameter  in  defining  the  response  of  a  structure  to 
high  intensity  acoustic  loading.  Increasing  this  parameter  provides  a  means  of 
reduoing  the  amplitudes  of  vibxatlon  and  associated  stresses  in  a  structure.  At 
present  no  reliable  method  for  predicting  damping  values  has  been  devised.  In  this 
Section  the  sources  of  damping  in  typical  structures  are  described  and  some  guidance 
on  the  range  of  values  likely  to  be  found  on  different  types  of  structure,  without 
special  damping  treatment,  is  provided.  Also  some  guidance  is  given  on  methods  of 
increasing  structural  damping. 

Damping  is  the  energy  dissipation  from  a  vibrating  system.  Two  theoretical 
models,  visoous  and  hysteretic,  are  used  to  represent  damping.  In  the  viscous  model 
the  energy  dissipated  per  cycle  is  proportional  to  the  product  of  the  frequency  and 
the  square  of  the  displacement  amplitude.  In  the  hysteretic  model  the  energy 
dissipated  during  a  harmonic  vibration  is  independent  of  frequency  and  proportional 
to  the  square  of  the  displacement.  For  structures  subjected  to  harmonic  excitation 
hysteretic  damping  may  be  treated  as  viscous  damping  using  an  equivalent  damping 
coefficient  inversely  proportional  to  the  driving  frequency.  In  this  Section  the 
viscous  damping  model  is  sed.  Many  terms  are  used  to  define  the  magnitude  of  damping. 
In  this  Section  the  measuxe  used  is  the  damping  ratio,  i.e.  the  ratio  of  the  actual 
damping  coefficient  to  the  critical  damping  coefficient  for  the  mode  of  vibration 
considered,  which  represents  the  total  damping  within  the  structure  from  all  sources. 
The  relationships  between  some  common  measures  of  damping  are  given  in  Appendix  2A. 

2.2  Damping  Sources 

2.2.1  General 


The  main  sources  of  damping  of  structural  response,  before  any  attempt  is 
made  to  Increase  it  by  artificial  means,  are  structural,  where  energy  is 
dissipated  within  the  structure  in  the  form  of  heat,  and  acoustic,  where  energy 
is  radiated  from  the  structure  as  sound.  Although  it  is  not  possible  to  give 
exact  relationships  between  the  magnitudes  of  structural  and  acoustic  damping 
ratios,  it  is  generally  expected  that,  for  typical  aircraft  structures, 
structural  damping  is  dominant.  Other  sources  are  present  but  generally  these 
make  relatively  small  contributions  to  the  overall  damping. 

2.2.2  Structural  damping 

Joints  and  material  internal  hysteresis  contribute  to  structural  damping. 
The  internal  damping  of  the  structural  material  is  very  low  in  conventional 
light  alloy  materials  used  in  aircraft  construction  compared  with  that  derived 
from  Joints. 

At  low  vibration  amplitudes  Joint  damping  results  from  hysteresis  in  the 
Joint  material.  As  the  vibration  amplitude  Increases  slip  occurs  in  the  area 
surrounding  the  rivets  and  damping  is  expected  to  increase  non-llnearly.  The 
onset  of  Joint  slip  depends  on  Joint  loading  and  the  normal  pressure  between  the 
Jointed  surfaces.  At  high  vibration  amplitudes  there  is  some  evidence  to  show 
that  Joint  damping  may,  Initially,  reduce  with  time  due  to  reduced  friction 
between  mating  surfaces. 

In  a  built-up  structure,  individual  joint  loadings  will  be  different  in 
each  mode;  hence  the  contribution  to  the  overall  damping  from  the  Joint  will  be 
different  in  each  mode. 

2.2.3  Acoustic  damping 

Both  overall  structural  vibration  modes  and  local  single  panel  modes 
contribute  to  acoustic  damping.  However,  it  is  not  possible  to  give 
quantitative  data  for  these  damping  values. 

Damping  from  modes  involving  distortion  of  the  whole  structure  is 
dependent  on  the  mode  of  vibration  and  the  wavelength  of  the  radiated  sound. 

Hie  acoustic  damping  ratio  for  a  typical  mode  varies  with  frequency  as  shown 
in  Sketch  (l).  When  the  wavelength  of  the  radiated  sound  is  greater  than  the 
longitudinal  wavelength  of  the  vibrating  mode  the  acoustic  damping  ratio  is 
saro. 
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Acoustic  damping  in  local  panel  modes  results  from  radiation  from  each 
element  of  the  panel  acting  as  a  piston  having  the  local  velocity  and  amplitude 
of  the  vibrating  panel.  When  adjacent  plates,  separated  by  a  stiff  edge  member, 
vibrate  in-phase  the  damping  pressure  on  one  plate  also  acts  on  adjacent  platrs 
increasing  the  overall  acoustic  damping  ratio.  If,  however,  adjacent  plates 
vibrate  in  anti-phase  their  damping  pressures  tend  to  cancel  each  other  giving  a 
small  overall  acoustic  damping  ratio.  The  variation  of  acoustic  damping  ratio 
with  plate  aspect  ratio  for  a  typical  rectangular  plate  vibrating  in  its 
fundamental  mode  is  shown  in  Sketch  (ll). 

2 . 3  Damping  Values  for  Typical  Structures 
2-3.1  Data  presented 

In  Figures  2. 1-2. 4  the  results  of  damping  measurements  made  on  different 
types  of  aircraft  structure  under  discrete  frequency  excitation  are  presented. 

In  each  figure  the  measured  damping  ratio  is  plotted  against  the  response 
frequency  without  differentiating  between  the  modes  of  vibration.  As  an  aid  to 
data  comparison,  where  sufficient  measurements  are  available,  a  line  is  drawn 
through  date  points  (using  a  least  squares  method)  assuming  that  the  damping 
ratio  is  inversely  proportional  to  frequency. 

2.3.2  Skin  and  stringer  panels 

In  Figures  2.1  and  2.2  measured  values  of  damping  ratio  are  presented  for 
flat  and  curved  panels  respectively.  These  figures  demonstrate  the  scatter  which 
is  typical  of  damping  measurements.  The  horizontal  shift  of  tha  mean  line  on 
Figure  2.2,  compared  with  Figure  2.1,  may  be  explained  by  the  increase  in 
frequency  due  to  the  curvature  of  th;  panels. 

The  damping  ratio  for  typical  aircraft  skin  and  stringer  panels  vibrating 
in  their  fundamental  mode  lies  in  the  range  0.005  to  0.030. 

2.3- 3  Integrally-machined  panels 

In  Figure  2.3  damping  ratio  is  plotted  against  frequency  for  built-up 
structures  having  integrally-machined  panels.  The  reduction  in  damping  ratio, 
compared  with  built-up  skin  and  stringer  panels,  is  attributed  to  there  being 
fewer  joints  in  machined  panel  structures. 

The  damping  ratio  for  a  typical  integrally-machined  panel  vibrating  in  its 
fundamental  mode  lies  in  the  ran,;*,  0.003  10  0.0i2. 

2.3- 4  Honeycomb  structures 

In  Figure  2.4  damping  ratio  is  plotted  against  frequency  for  flat 
honeycomb  panels  and  honeycomb  wedges 

The  damping  ratio  for  a  typical  hone>  .-omb  panel  vibrating  in  its 
fundamental  mode  lies  in  the  range  0.010  to  0.030. 


Factors  Affecting  Damping 
2.4-t  Structural  .joint* 

The  area  of  contact  between  jointed  surface*  and  the  noreal  pressure 
between  them  affects  the  Joint  load  at  which  slip  first  occurs;  damping  is 
therefore  dependent  on  the  joint  efficiency. 

The  use  of  a  viscoelastic  Jointing  compound  can  modify  the  atresa 
distribution  near  the  joint  and  reduce  fretting  but  might  not  increase  the 
overall  damping  ratio  (see  Figure  2.1).  The  separation  of  the  metal  surfaces 
by  the  Jointing  compound  tends  to  prevent  any  reduction  in  damping  with  time 
resulting  from  oxidation  of  rubbing  metal  surfaces. 


2.4.2  Frequency 

Viscous  damping  ratios  are  inversely  proportional  to  the  frequency  of 
vibration  and  hysteretlc  damping  ratios  are  inversely  proportional  to  the  square 
of  the  frequency  of  vibration.  As  each  individual  damping  source  is  most 
closely  represented  by  one  of  the  two  damping  models,  an  assemblage  of  all 
damping  sources  usually  involves  a  combination  of  both  models  so  that  a  general 
relationship  between  damping  and  frequency  is  not  possible. 

2.4.3  Skin  Thickness 

The  effect  of  skin  thickness  on  the  damping  ratio  is  dependent  on  the 
mode  of  vibration.  For  a  stiffened  panel  vibrating  in  its  stringer  bending  mode, 
damping  increases  with  plate  thickness;  as  plate  thickness  is  increased  damping 
in  the  stringer  torsion  mode  passes  through  a  maximum.  Further  information  on 
the  effect  of  skin  thickness  on  damping  is  given  in  Derivation  2. 7. 4* 

2.4.4  Curvature 

Limited  experimental  data,  on  similar  structural:,  suggest  that  curvature 
has  little  effect  on  the  damping  ratio.  Curvature  is  likely  to  have  the 
greatest  effect  in  panel  modes  where  stringers  bend.  An  approximate  method  for 
estimating  the  effect  of  curvature  on  the  damping  ratio  for  a  particular  mode  ie 
given  in  Derivation  2-7.4. 

2.4.5  Flow  velocity 

Air  flowing  past  panels  can  reduce  the  damping  ratio.  Thi#  ie  important 
in  aircraft  structures  where  the  reduction  in  damping  may  lead  to  panel  flutter 
particularly  in  supersonic  and  transonic  flow  regions . 

2.4.6  Amplitude  of  vibration 

The  variation  of  damping  with  amplitude  of  vibration  for  a  typical 
excitation  is  shown  in  Figure  2.5.  All  data  presented  on  this  Figure  are  from 
one  series  of  tests  in  a  travelling  wave  tube.  Damping  ratios  calculated  using 
the  single  mode  (all  edges  fixed)  response  theory  are  plotted  against  the  root 
mean  square  panel  st-sss.  Thors  is  some  indication  that  damping  ratios  are 
reduced  with  increase  in  response  amplitude  for  these  specimens,  which  la  against 
the  expected  treno  for  structural  damping  (see  Section  2.2.2).  In  Reference 
2.7.13  data  for  aar.dj'ich  panels  show  an  increase  In  damping  ratio  with  increase 
in  stress  amplitude. 

Methods  of  Increasing  Damp in 


In  general,  daaiping  treatments  have  to  far  only  been  used  on  a.’rcraft  structures 
to  tllevlate  a  service  problem  not  foreseen  at  the  design  stage.  A  better 
undt rs tending  of  the  advantages  of  good  damping  should  eventually  lead  to  its 
consideration  as  a  design  criterion.  Consideration  of  damping  at  the  design  stage  may 
allow  panels  with  sandwiched  viscoelastic  material  to  be  used  which  era  continuous 
across  atrlnger  joints.  However,  where  primary  structure  it.  concerned,  there  are 
difficulties  to  be  overcome  including  loss  of  stiffness  in  Joint  fabrication  and  creep 
under  rteady  load. 

To  increase  damping  in  a  structure  efficiently,  a  knowledge  of  the  frequency  end 
shape  of  the  modes  requiring  additional  damping  is  necessary.  Changing  law  stiffness 
of  a  structure  does  not  improve  damping  unless  the  critical  resonant  frequency  can  be 
made  to  approach  the  frequency  for  peak  acoustic  damping.  However  increasing 
structural  stiffness  Increases  the  natural  frequencies  and  may  thereby  provide  an 
acceptable  remedy  for  excessive  vibration  if  the  higher  frequencies  are  not  sxclted 
during  normal  oporatlng  conditions. 


H 


The  us*  of  viscoelcetic  damping  material  provides  the  Moat  practical  method  for 
improving  damping.  To  be  effective  these  materials  must  bo  applied  so  thr t  they  will 
undergo  high  strain.  It  is  alt’o  necessary  to  select  a  material  with  suitable  physical 
properties  ovar  the  range  of  operating  temperatures.  hat* rial  properties  vary  great!/ 
with  temperature  and  are  also  frequency  dependent.  Th<  variation  of  material 
properties  with  temperature  for  a  typical  viscoelastic  daa<ping  compound,  intsndsd  for 
use  at  normal  ambient  temperatures,  is  shown  in  Skstcb  (ili). 


TW»p«roti,'r*  °C 


Viscoelastic  material  may  be  applied  to  a  structure  as  an  unconstrained  layer. 

For  optimum  damping  effect  the  material  should  be  placed  in  a  region  of  *  gb  bending 
strain.*  The  regions  of  highest  bending  strain  are  at  the  anti-node  positions,  and  it 
is  an  advantage  to  have  the  damping  material  as  far  removed  as  possible  from  the 
neutral  axle. 

Unconstrained  viscoelastic  material  need  only  be  used  In  areas  of  high  bending 
strain  in  the  modes  requiring  additional  damping.  However,  if  several  modes  need 
increased  damping  it  may  be  necessary  to  cover  the  whole  skin  surface.  For  stiffened 
panels  having  stiffenera  jf  low  flexural  stiffness,  such  as  integrally-machined  panels, 
it  is  an  advantage  to  use  viscoelastic  material  on  the  stiffener  face  remote  from  the 
skin. 


Constrained  viscoelastic  layers,  formed  by  sandwiching  damping  material  between 
skins,  give  increases  in  damping  when  the  viscoelastic  material  ie  subjected  to  high 
shear  strain.  A  similar  constrained  layer  damping  system  is  formed  by  attaching  foil- 
backed  damping  tape  tc  the  skin  surface.  Expression •  and  calculated  values  for 
damping  ratios  of  panels  with  viscoelastic  damping  layers  are  given  in  Derivation  2.7-7 
and  Reference  2.7.14. 

Damping  in  structures  may  be  increased  by  using  viscoelastic  inserts  at  joints. 
This  method  is  not  likely  to  give  as  high  an  increase  in  damping  ratio  as  plate-surface 
layer  treatments  but  it  reduces  joint  fretting. 

Addition  of  viscoelastic  material  leads  to  an  increase  in  mass  and  changes  in 
stif^neas  the  structure.  These  changes  should  be  considered  in  optimising  the 
overall  effect  of  the  damping  materials. 

2.6  Measurement  of  Damping 

Damping  is  commonly  measured  fro«-  resonant  response  peak*,  logarithmic  decrement 
or  from  vector  plots  (see  Reference  2. 7.12).  Vhen  natural  frequencies  are  not  well 
separated  care  is  necessary  to  ensure  thet  ths  measured  damping  Is  that  for  the  mode 
required  and  not  a  combined  damping  ratio  from  a  number  of  modes.  To  measure  damping, 
instruments  such  as  strain  gauges  are  used  which  are  connected  by  cables  to  recording 
equipment.  The  instrumentation  and  cables  add  to  the  mass  and  stiffness  of  the 
structure  and  may  also  cause  artificial  increases  in  damping.  Because  of  these  and 
other  experimental  difficulties,  measured  damping  values  should  be  considered  as 
upper  boundary  values |  this  is  particularly  so  when  damping  is  measured  in  siren 
facilities  when  energy  dissipation  along  the  duct  is  likely  to  result  in  r.n  overestimate 
of  the  true  damping. 

2 . 7  Derivation  and  References 
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FIGURE  2  3.  8UILT-UP  STRUCTURES  WITH  INTEGRALLY  MACHINED  SKINS 
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FIGURE  2.5.  AMPUTUOE  EFFECT  ON  DAMPING  RATIO 
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APPENDIX  2k 


2A. 1  Notation 


lbf  s/ in 


lbf  s/in 


lbf/ in 
lbf/ in 
lbf  a^/in 


l  amplitude  constant  * 

3  damping  coefficient  N  s/m 

3  oritical  damping  coeff  loieut,  i.e.  value 

°  of  damping  coefficient  for  which  a  fre» 

motion  is  just  aperiodic  M  s/m 

if  width  of  frequency  response  curve  at  which 

amplitude  is  1//2  times  maximum  amplitude  Hz 

P(t)  time  dependent  force  N 

Po  amplitude  of  forcing  function  N 

f  frequency  of  vibration  Hz 

rn  natural  frequency  of  vibration  Hz 

tt  damping  component  of  hysteretic  stiffness  N/m 

Ic  spring  stiffness  N/m 

a  mass  of  vibrating  system  leg 

3  magnification  factor 


dlaplacamant 


A  logarithmic  decrement 

6  damping  ratio 

6jj  hysteretic  damping  ratio 

ij  loas  factor 

2A.2  Damping  Ratio 

Considering  the  single  degree  of  freedom  system  with  a  viscous  damper, 
represented  dlagraamatically  in  Sketch  2A(i),  the  equation  of  motion  of  the  free  end  of 
the  spring  relative  to  the  unatrained  position  is 

m  +  o  —  +  kx  -  P(t)  (l) 

dt^  dt 


Sketch  2A(1) 

The  solution  for  free  vibration  (when  F(t)  *  0)  for  a  lightly  damped  system 
( 6<1 . 0)  is 

x  *  Aq  cos(2*fnt y  l-ft*) ,  (2) 

where  6  is  the  damping  ratio,  •nd  °crit  ls  givan  by  2,/mk. 
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2 A. 3  Magnification  Factor 

When  r(t)  represents  the  harmonically  varying  force  F^cos2  ft,  the  ratio 
| x j /F  is  oqs&l  to 


Plotting  |x|  /F  against  f  gives  the  frequency  response  curve  shown  In  Sketch  2A(il). 


The  magnification  factor,  Q,  is  the  ratio  of  the  maximum  displacement  (resonance 
value)  to  the  static  force  (f=o)  displacement.  From- Expression  (3)  it  is  evident 
that  q  »  1/26. 

Provided  that  6  is  not  greater  than  about  0.3  it  car  be  shown  that 

-  26  =  A  . 

fn  “ 


2 A. 4  Logarithmic  Decrement 

/ 

The  period  of  the  damped  free  oscillation  of  Expression  (2)  is  2 x/  2*f«/l-6  ; 

n/ - j 

_ «  f  /  — 

hence  the  ratio  of  maximum  amplitudes  of  successive  cycles  is  given  by  e  '  * 


The  logarithmic  decrement, A  ,  is  given  by  the  natural  logarithm  of  the  ratio  of 
two  successive  amplitudes.  The  amplitudes  measured  at  the  point  of  tangency  betwoen 
the  oscillation  envelope  and  the  decaying  oscillation,  shown  In  Sketch  2A(ijl),  closely 
approximate  to  the  ratio  of  maximum  amplitudes  of  successive  cycles.  For  lightly 

damped  systemo  where  J 1  -  62  approximates  to  unity  A  =  2*5. 


kotch  2A( ill) 


2A 


Loss  Factor 


Th*  previously  defined  duping  values  are  based  on  a  viscous  damping  system. 

For  harmonically  excited  systems  a  complex  stiffness  type  of  damping  can  be  used.  The 
stiffness  and  damping  forces  are  represented  by  (k  +  ih)x.  This  expression  can  be 
written  is  (it  1 7 )  lot  where  7  is  the  less  factor  (7  *  h/k)  . 

For  this  type  of  damping 


Putting  «  h/2k  this  expression  reduces  to  the  viscous  damping  form. 

2 A. 6  Damping  Measure  Relationship 

Th  1  relationship  between  the  viscous  damping  ratio  and  other  damping  measures 
for  harmonic  vibration  is  summarised  in  the  table  below. 


Viscous  dan  ping 
ratio 


Magnification 

factor 


Logarithmic 

decrement 


Loss 


factor 


6 


_1_ 

2Q 


_A_ 

2* 
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Section  3 

REFERENCE  FREQUENCY  OF  PANEL  WITH  FLEXIBLE  STIFFENERS 


3. 1  Notation 


A. 

cross-sectional  area  of  stiffener  (see  Sketch  (i)) 

2 

m 

in2 

A 

length  of  longer  side  of  plate 

m 

in 

b 

length  of  shorter  side  of  plate 

m 

in 

D 

flexural  rigidity  of  plate  given  by  Ept3/{l2(l->o-2)j 

N  m 

lbf  in 

E 

P 

Young's  modulus  of  plate  material 

N/m2 

lbf/in‘ 

E. 

Yo'ing's  modulus  of  stiffener  material 

N/m2 

lbf/ in‘ 

f 

reference  frequency  of  plate 

H* 

c/s 

I 

second  moment  of  area  of  stiffener  (without  skin) 
aoout  its  neutral  axis  parallel  to  plate 
(see  Sketch  (i)) 

«* 

in1* 

K 

reference  frequency  parameter 

m/e 

in/  s 

t 

thic’uib's  of  plate 

m 

in 

V 

velocity  parameter  for  plate  material” 

or 

Poisson's  ratio  for  plate  material 

P 

density  of  plate  material 

kg/»3 

♦ 

P* 

density  of  stiffener  material 

kg/m3 

Both 

SI  and  British  units  are  quoted  but  any  coherent  system 

of  units 

may  be 

used* 

*  The  velocity  parameter  ia  defined  in  SI  unite  by 

V  =  (E/p)1/?/5080  or  in  British  unite  by 

V  *  (E/p)1,/2/200  OOO.  It  io  approximately  unity  for  all  common  structural 
metallic  materials. 

+  A  density  value  expressed  in  British  units  as  pounds  per  cubic  inch  has  to  be 
divided  by  3 86.4  before  being  used  to  calculate  parameters  defined  in  this  Section. 

2 

(A  forco  of  1  lbf  acting  on  a  mass  of  1  lb  produces  an  acceleration  of  386./*  in/s  .) 


S*tin 


A*i>  lor  I 


Stiftow 


At  0 ed  I  rslota  to 
•haded  oreo  only 


Sketch  (1) 


3.2  Notes 

This  Section  gives  the  reference  frequency  of  initially  unstressed,  flat,  uniform, 
rectangular  plates  in  a  panel  with  flexible  stiffeners.  The  reference  frequency  is 
the  fundamental  natural  frequency,  obtained  from  the  model  used  to  represent  the  panel, 
where  an  assumed  nude  shape  has  been  employed  in  a  Tayleigh-Ritz  analysis.  This 
frequency  may  be  used  as  a  first  approximation  to  the  fundamental  natural  frequency  of 
the  plate  with  edgei  fixed  against  rotation.  The  reference  frequency  is  intended  for 
use  in  calculating  the  stress  response  of  panels  with  flexible  stiffeners  under  the 
action  of  acoustic  loading.  In  the  vibration  mode  considered  the  stiffeners  nre 
assumed  to  be  free  to  bend  in  a  plane  normal  to  the  plate  and  be  fixed  against  rotation 
I  (see  Sk-'-ch  (  ii )  )  . 

1 

I 

! 

i 
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The  reference  frequency  ia  given  by 


In  Figures  3*l-3*4  value*  of  K  are  plottad  against  E#l/Db  for  a  range  of 
values  of  a/b.  Each  of  these  figures  is  drawn  for  a  particular  value  of  p^A^/pbt. 

The  plates  are  assumed  to  be  separated  by  flexible  stiffeners  along  their  longer 
sides  and  fully  fixed  along  their  shorter  sides.  In  the  derivation  stiffener  bending 
is  considered  but  stiffener  toreion  is  neglected. 

The  use  of  a  second  moment  of  area  of  the  stiffener  (without  skin)  about  its 
neutral  axis  parallel  to  the  plate  has  been  found  to  give  closest  correlation  with 
limited  test  data.  It  is  possible  with  cex'tain  configurations  that  the  effective 
neutral  axis  could  lie  closer  to  the  skin  but  as  yet  there  is  insufficient  evidence 
to  clarify  this  point. 

This  Section  is  particularly  suitable  for  the  estimation  of  the  natural  frequency 
of  integrally-machined  panels. 

In  deriving  the  curves  the  value  of  <r  was  assumed  to  be  0.3*  This  value  gives 
sufficiently  accurate  frequencies  for  all  common  structural  metallic  materials. 

A  computer  program,  based  on  this  procedure,  to  calculate  both  the  reference 
frequency  and  stress  response  is  given  in  Section  4  of  this  AGARDograph. 

3*3  Comparison  with  Measured  Data 

Figure  3.5  shows  a  comparison  of  calculated  reference  frequency  and  lowest 
measured  natural  frequencies  of  flat  stiffener  panels.  For  some  of  the  integrally- 
machined  panel  frequency  data,  where  measured  frequencies  are  taken  from  strain 
response  test  data,  frequency  bands  covering  the  predominant  response  frequencies  atre 
shown. 


It  should  be  noted  that  in  estimating  frequency  the  plate  edge-fixing  approximates 
to  the  fully-fixed  condition  since  all  edges  are  fixed  against  rotation  and  only  the 
longer  plate  edges  are  free  in  translation.  Vith  the  restriction  ou  plate  edge 
rotation  the  estimated  frequencies  are  expected  to  be  greater  than  the  lowest  measured 
values  which  correspond  with  modes  involving  stringer  torsion. 

3*1*  Derivation 
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3«5 

It  la  required  to  estimate  the  fundamental  natural  frequency  of  a  flat  aluminium 
alloy  plate  in  a  panel  array*  The  plate  hat  the  following  dimensions  and  physical 
properties  t 

A#  -  43*3  mm2,  I  «  950  mm4, 

a  ■  600  mm,  t  >0*8  mm, 

b  >  120  sub,  7  ■  0*3, 

£  «  Eo  -  70  000  MN/m2,  P  *  P.  '  2770  kg/m3. 


Firstly 


70  000  x  106  x  (0.8  x  IQ'3)3 

12(1  -  0.32) 


-  3*28  N  m, 


2770  x  43.3  z  10 
2770  x  120  x  0.8  x  10" 


-  0.451, 


70  000  x  10 


50  x  10 


3.28  x  120  x  10"3 


From  Figure  1  for  paAa/pbt  =0  K  *  2680  m/s 

from  Figure  2  for  p^A^/pbt  =1.0  K  *  1980  m/s 

and  from  Figure  3  for  paAg/pbt  *  2.0  K  *  I67O  m/s 

therefore  by  interpolation  for  PgA^/pbt  *  0.451  K  =  2260  m/s. 


from  Figure  2  for 
and  from  Figure  3  for 


70  000  x  10 


Hence  the  fundamental  natural  frequency  of  the  panel  is 


0.99  x  2260  x 


0.8  x  10 
1202  x  10" 


124  Hz 


r«re  pr^'Y*'  r'wv*.  «%■»•*»*»  — -  «*•»  »mwrrtwefi>t,y  ws-  8,w*,>?«p»w<  Hf  <W(W!^*W5 


FIGURE  3.3.  FREQUENCY  PARAMETER  FOR 


r\-ry  “rrjj'  ' 
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Section  k 

THE  ESTIMATION  GF  R.M.S  STRESS  IN  SKIN  PANELS  WITH  FLEXIBLE 
STIFFENERS  SUBJECTED  TO  RANDOM  ACOUSTIC  LOADING 


4.1  Notation 


A. 

cro*i-iec tional  ».'»»  of  stiffener 
(•••  Sketch  (i)) 

«2 

in2 

a 

.length  of  longer  side  of  plate 

m 

in 

b 

length  of  shorter  aide  of  plat* 

■ 

In 

D 

plate  flexural  rigidity  given  by 

EptVflZd-o2)} 

N  m 

lbf  in 

d 

height  of  etlffener 

m 

in 

E 

P 

Young's  modulus  of  plate  material 

»/■* 

lbf/ in2 

£ 

e 

Young's  aodulus  of  stiffener 
material 

N/m2 

lbf/ in2 

f 

fundamental  natural  frequency  of 
panel  with  edges  fixed  againat 
rotation 

Hz 

c/s 

Op(f) 

spectral  density  of  acoustic 
pressures  at  frequency  f 

(N/il2)2/H* 

(lbf/in2)2/ (c/s) 

i 

second  moment  of  area  of  stiffener 
(without  akin)  about  its  neutral 
axis  parallel  to  plate  (aee  Sketoh  (i)) 

X* 

stress  paramt  ,er  (for  plate  atreeeea 

K'*K  and  for  stiffener  stress 

n 

KME^/EpMd/tJX^) 

X 

n 

plate  stress  parameter  (n«l- 5  and 
denotes  location  of  plate  stress  - 
see  Sketch  (ii)) 

K 

• 

stiffener  stress  parameter 

K6 

damping  ratio  correction  factor 

Lp-(f) 

spectrum  level  of  acoustic  pressure 
at  froquency  f 

dB 

dB 

P 

uniform  static  pressure  on  plate 

N/m2 

lbf/in2 

p 

*rms 

root  mean  square  fluctuating 
pressure 

N/m2 

lbf/in2 

S 

0 

ratio  of  stress  at  considered  plate 
location  to  applied  uniform  static 
pressure  on  plate 

S 

rms 

r.m.Sr  stress  at  sviffener  due  to 
acoustic  loading 

N/*2 

lbf/in2 

t 

thickness  of  plate 

■ 

in 

6 

damping  ratio  in  vibrating  mode 

6ref 

referenoe  damping  ratio  (*0.017) 

3i 


p  density  of  plats  Material 

Pa  density  of  stiffener  Material 

0  Poisson's  ratio  for  plats  Material 


kg/m3  ♦ 

kt/M3  + 


Stan 


A,  ond  I  ratal*  to 
SwSad  oraa  only 


SMch  (i) 


Both  SI  and  British  units  are  quoted  hut  any  oohsrsnt  systSM  of  units  My  be  used. 

+A  density  value  expressed  in  British  unite  as  pounds  per  oubic  inch  has  to  be 
divided  by  386.4  before  being  used  to  calculate  paraaeters  defined  in  this  Section. 

(A  force  of  1  lbf  acting  on  a  Mass  of  1  lb  produces  an  acceleration  of 

386.4  in/*2.) 

4.2  Notes 

This  Section  gives  a  Method  of  estimating  the  r.u.s.  stress  in  flat  rectangular 
skin  panels  with  flexible  stiffeners  under  the  aotlon  of  randoM  acoustic  loading. 

The  stress  locations  on  the  plate  and  stiffener,  and  directions  of  r.m.s.  stress, 
for  which  values  are  provided  are  shown  in  Sketch  (ii). 


The  r.a.s.  stress  for  a  stiffened  panel  subjected  to  randoM  acoustic  loading  on* 
side  is  given  approximately  by  the  expression 

1/2 

4.1 


ras 


'ft 

ref  y 


X.  S  , 
o  o  ’ 


when  only  the  fundaaental  natural  mode  is  excited  by  the  noise.  Por  plate  etressoe 
the  ratio  SQ  is  proportional  to  the  product  of  a  stress  parameter  and  (a/t)  and 
for  the  stiffener  stress  SQ  is  proportional  to  the  product  of  a  stress  parameter, 

O 

(a/t)  and  (d/t).  In  Figures  4.1-4.20  values  of  KQ  are  plotted  against  I^I/Db 
for  a  range  of  values  of  a/b.  Similarly  in  Figures  4.21-4.24  values  of  JCa  are 
presented.  Bach  of  these  Figures  is  drawn  for  a  particular  value  of  paAa/pbt  as 


detailed  in  Table  4.1 


Tk’»n«T?*> 


■r  *r ■vrrjvir  "“ 
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TABLE  4,1 


Figure  nos.  for  stress  paramstsr 

Vs 

Strese^-^ftot 

location 

0 

1 

2 

3 

1 

4.1 

4.2 

4.3 

4.4 

2 

4.5 

4 .6 

4.7 

4.8 

3 

4.9 

4.10 

4.11 

4.12 

4 

4.13 

4.14 

4.15 

4.16 

5 

4.17 

4.18 

4.19 

4.20 

S 

4.21 

4.22 

4.23 

4.24 

Figures  4.25  end  4.26  give  nomographs  for  srMS*  Th»  nomograph#  are  entered  at 

a  value  of  Lpa(f),  each  quadrant  being  used  in  turn  in  the  direction  indicated  through 

ranges  of  a/t,  K'  and  f.  Figure  4. 26  is  an  extension  of  the  stress  range  of 
Figure  4.25.  The  stress  nomographs  “.re  drawn  for  a  value  of  6  a.  '>.017}  for 

other  values  of  b  is  calculated  using  a  correction  factor  whiv \  is  plotted 

against  b  in  Figure  4.27.  Guidance  on  values  of  b  for  typical  at  -aft  structures 
say  be  found  in  Section  2  of  this  AGARDograph. 

The  use  of  a  second  moment  of  area  of  the  stiffener  (without  ••icln)  abi  t  its 
neutral  axis  parallel  to  the  plate  has  been  found  to  give  closest  .'rrela tio.<  with 
limited  test  data.  It  is  possible  with'  certain  configurations  that  the  effective 
neutral  axis  could  lie  closer  to  the  skin  but  as  yet  there  is  insufficient  evidence 
to  clarify  this  point. 

For  panels  with  riveted  stiffeners  the  edge  stresses  may  be  used  as  the  nominal 
unfactored  rivet  line  values.  In  using  these  stresses  to  calculate  a  fatigue  life 
care  must  be  taken  to  ensure  that  fatigue  data  from  a  similar  structure  aro  used  so 
that  geometric  stress  concentration  factors  are  of  the  same  order. 

In  obtaining  values  of  K  and  Kg  presented  in  Figures  4,1-4.24  the  panels 

were  assumed  to  be  built-up  from  plates  fully  fixed  on  their  shorter  sides,  and 
attached  to  flexible  stiffeners  along  their  longer  sides.  The  assumed  deflected 
shape  is  shown  in  Sketch  (ill).  At  plate  locations  1  and  2  the  value  of  Kn ,  for 

plates  having  high  aspect  ratio  and  very  flexible  stiffeners,  reduces  to  zero  and 
then  increases  as  E^I/Db  increases.  This  is  due  to  the  change  in  sign  of  tne 

static  stress  ratio  as  the  relative  stiffness  of  the  plate  to  that  of  the  stiffener 
increases.  The  number  of  individual  plates  in  the  panel  array,  between  each  pair  of 
rigid  supports  or  stiff  frames,  was  taken  to  be  sufficiently  great  for  the  centre 
panel  of  the  array  to  be  unaffected  by  tne  number  of  plates  In  the  panel. 


Skatch  (hi) 
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The  pressure  ia  assumed  to  bo  uniform  and  in  phaaa  ovar  tha  whole  of  each 
individual  plate  and  tha  spectrum  level  of  tha  aooustio  pressure  ia  taken  to  ba 
conatant  ovar  tha  range  of  fraquenclaa  cloaa  to  the  fundamental  natural  frequency 
of  the  panel. 

In  tha  derivation  of  the  nomograph  it  ia  aaaumed  that  the  plate  banding  atreaa 
ia  within  tha  linear  region  where  it  ia  directly  proportional  to  the  normal  pressure, 

that  ia  p/Ep  leaa  than  about  20(t/b)^  (aee  Reference  4.3.10).  Also  the  value 

of  6  ia  taken  tc  be  0.3;  uae  of  thia  value  gives  a  sufficiently  accurate  stress 
for  all  common  structural  metallic  materials. 

For  conventional  structures  without  special  damping  treatment  f  should  be 
taken  aa  the  undamped  natural  frequency. 

This  Section  ia  particularly  applicable  to  integrally-machine  panel  structures. 

In  Appendix  4A  a  computer  program,  based  on  the  procedure  of  this  Section  and 
Section  3,  ia  described  which  calculates  the  fundamental  natural  frequency  and  stress 
response  to  acoustic  loading  of  panels  with  flexible  stiffeners.  Sub-routinas  are 
given  for  both  the  frequency  and  stress  calculations. 

4*3  Calculation  Procedure 

4.3.1  The  procedure  t-  estimating  S  in  the  general  case  is  as  follows. 

nn» 

(i)  Estimate  the  reference  frequency  of  the  panel  using  Section  3  of  this 
AGARDograph,  This  reference  frequency  may  be  taken  as  the  fundamental 
natural  frequency  for  the  purpose  of  calculating  stress. 

(ii)  Obtain  the  value  of  spectrum  level  of  acoustic  pressure  Lpa(f)  at  the 

calculated  frequency.  If  only  the  band  pressure  level  is  known,  it  is 
first  corrected  to  pressure  spectrum  level  (unit  bandwidth)  using 
Reference  4.3.8.  The  reference  pressure  for  sound  pressure  level  is 

20  pN/m2. 

(iii)  Calculate  the  parameters  a/b,  EaI/Db  and  p^A^/pbt  and  from  Figures 
4.1-4.24  obtain  the  approprate  value  of  K'  or  Ka. 

(iv)  For  stiffener  S  calculate  the  parameter  (E  /E  )(d/t)  and  factor 

ms  s  p 

by  that  parameter  to  obtain  the  required  value  of  K*  . 

(v)  Calculate  the  parameter  a/ 1  and  read  the  value  of  Srja8  from  Figure 

4.23  or  4.26.  The  nomographs  ars  entered  at  a  value  of  L  (f),  each 

ps 

quadrant  being  used  in  turn  in  the  direction  indicated  through  ranges  of 
a/ 1 ,  K*  and  f. 

(vi)  For  values  of  6  other  than  0,017,  factor  the  estimated  value  of  Srjna 
by  Kq  obtained  from  Figure  4.27. 

4.3.2  Within  the  nomograph  the  spectrum  sound  pressure  level  is  converted  into 
the  spectral  density  of  acoustic  pressure.  The  spectrum  sound  pressure  level 

2 

is  converted  into  the  root  mean  square  fluctuating  pressure  in  units  of  (N/m  )/Hz 
(tee  expression  below  or  Reference  4.3.9)  and  then  squared  giving  a  value  in 

2  2  2 

units  of  (N/m  )  /Hz  ,  Since  unit  bandwidth  ia  used  this  ie  numerically  equal 

2  2 

to  the  spectral  density  of  acouetic  pressure  Gp(f)  in  units  of  (N/m  )  /Hz. 
Lp.(f)  «  20(log10  prs!9  +  4.70). 

If  (f)  is  required  in  British  units  of  (lbf/in2 )2/(o/s)  it  ie 

given  by  8 

Lp.(f)  -  2°(l°*10  Prms  +  8*54)’ 
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4.4  Comparison  with  Measured  Data 

Figure  4*28  show*  a  comparison  of  estimated  and  measured  stresses  xn  flat 
stiffened  panels  using  the  calculated  frequency  and  assuming  6  *  0.017, 

Figure  4.29  shows  a  comparison  of  estimatod  and  measurrd  stresses  using  average 
measured  values  for  f  and  6. 
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4.6  Example 

It  is  required  to  estimate  the  plate  surface  stresses  at  the  centre  and  eiges 
on  the  plate  centre  line,  and  the  stiffener  tree-edge  stress  at  the  centre  of  the 
longer  plate  side.  The  plate  is  an  element  of  an  internally-machined  panel 
subjected  to  Jet  noise  on  one  side. 

The  variation  ox'  sound  pressure  level  over  a  range  of  frequencies  is  given  in 
the  table,  sound  pressure  levels  being  1/3  octave  band  levels. 


Sound  pressure 
,  ,  /  \ 
xevcj.  vaiJ/ 

I 

151 

i  55 

156  1 

' 

135 

153 

Frequency  (Hz) 

50 

100 

200 

300 

- - - 

500 
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The  panel  is  made  up  from  uniform  plates  having  the  following  dimensions  and 
properties : 


b 

d 

E 


43.3  nsn2, 

600  mm , 
120  mm, 
16.2  mm, 


E 


70  000  MN/m  , 


I 

t 

0 

& 

p 


950  nmH, 
0.8  mm, 

0.3, 

0.015, 


*  2770  kg/ nr  . 


From  Section  3  the  fundamental  natural  frequency  of  the  panel  is  124  Has. 

By  interpolation  from  the  table  the  1/3  octave  band  pressure  level  at  124  Hz 
is  155.9  dB. 


From  Reference  4.5-8, 

L  (f)  =»  155.9  -  14.8 

po 


141.1  dB. 


Kow 


u 

=3 

12  <1-0. 32) 

a  j  *  to  1 

a 

b 

= 

600 

*  s 

120 

5. 

E  I 

a 

70  000 

x  106  x  950  x  10‘ 

-12 

=  169 

Db 

3. 

28  x  120  x  10~3 

P  A 

a  a 

2'. '70 

x  43.3  x  10“6 

=  0.451 

pbt 

2770  X 

120  x  0.8  x  iO"6 

a 

t 

a 

600 

0.8 

7  50  and  —  s 
t 

16.2 

0.8 

=  20.2. 

and 


Also 


From  Figure  4.27,  for  6  =  0,015,  =  1.06. 

Plate  r.m.s,  stresses 

From  Figures  4.1-4.20  the  values  of  K'  are  read  for  p^A^/pbt  =  0  and  1  and 

by  interpolation  the  values  of  K'  for  PgAB  =  0.451  are  found.  From  Figure 

4.26,  entering  the  nomograph  at  141.1  dB  the  values  of  plate  r.m.s.  stresses  are 
obtained  as  shown  in  the  table. 


Stress  parameter  K' 

HD 

Panel 

location 

P  A 

■  '■*  «  0 

Pbt 

^.1 

fbt 

P  a 

—  -  0.451 
pbt 

|  ~r  mam  |  hhmi 

H 

D 

1 

.0264 

.0135 

,0206 

121 

128 

1  2 

.0053 

.0009 

0033 

19.4 

20.5 

3 

.0040 

.0037 

.0039 

22.9 

1 

24.3 

4 

.0122 

.0080 

.0103 

6  0.5 

64.I 

5 

.0285 

.0156 

.0227 

i33 

141 
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Stiffener  |tm» 

From  Figures  4. 21-4.24  the  value  of  K  is  found  as  for  the  plate  stresses. 

Factoring  by  (E  /E  /( d/t)  gives  the  required  value  of  K.' . 
s  p 

The  value  of  r.ni.s.  stress  is  found  from  Figure  4.26. 


Stress  parameter 

E  d 

K'  .  K  7 

s  E  t 

P 

Sr»s  MN/"2 

P  A 

o  o 

Pbt 

P  A 

-a-a  -  i 

pbt 

P  A 

-2-i  =  0.451 

Pbt 

6  c  0.017 

6  «  0.015 

.0027 

.0028 

.0027 

.0562 

330 

350 

I 
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FIGURE  4.15.  STRESS  PARAMETER  AT  LOCATION  4  FOR  />,  A, //>„ bt  =  2 


k6 


hi 


10  \d  e,i  io5  io4 

w 

FIGURE  4.21.  STRESS  PARAMETER  AT  LOCATION  S  FOR  p1At/^>pbt  =  0 


Ob 

FIGURE  4.22.  STRESS  PARAMETER  AT  LOCATION  S  FOR  p,A,//>pbt  =  l 


■MiaMiai 
■  ■■■■■■I 

•■■Mini 


n/i'MM'j 


MnUHIlIHr/WKflrillll 

ggfjjaBBMIl  I 


■■■■■■■■ 

■■■mi 

■■11111 
[Hill 


■  WSJ  »* 

EBmuna^am 

-  .  <111111 
^nBiEiiiii 


aGBllll 

■limn 

mm 

■1 

I 

1 

SSSSi'iiiSnS 

— ■■■■MIHMM 

mm 

E 

*■1 

K 

iiM'arii 

■'ilMI 

nuiii 

■bihi 

Esauu 

■■nun 
liiiiiii 
■■■Hill 
111111111 


■■■■■■■a 

hmiiiih 

lAllilll 

■■■■llll 

lillllll 

■min 


■  aaania 


MBsaMflaaaal 

IB  ■■■■■■■■ 
■■■■■■■III 
■■■■■■■11 

M 

■  manaBail 

■■■■■■■III 

■■■■■lii!| 

1111111 


|5mMbnr^riiiMMMiKPiniigHaaaaanMi»nrinMiniai{ 

l<&Kmmimn%t%n’mzEtKmarJitwimm’BtmamntMmrAWJw,Mmnu\ 

\wmrMwm^^mwsmnnwmmummwMMmrAmw:i\\i 

mmmm&im  mmwmsmm 

bbmhbIhmimmmmmm 


I  mmVM 

a.ar.ai’.M — — ^ ^ 

■bMKa'ii»iim!iMagtiiiir«Riuu^ariraMM  ■■■■»» 

■■■rjaiiuMUiHiiiinsioiiaNpaHiiHiii 

■■■Miiiii  nBBianiiHBaiiiii 


ini  mmamimmm 

lililaiiiiiHiiiii 


ii’ViitiiUi 


Estimated  S— -  MN/m 


52 


f 

i 


FIGURE  4.28  COMPARISON  OF  ESf (MATED  ANO  MEASURED  STRESS  USNG  CALCULATED 
FREQUENCY  AND  8  =  0  017  FOR  STRESS  ESTIMATION 
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FIGURE  4.29  COMPARISON  OF  ESTIMATED  AND  MEASURED  STRESS  USING  MEASURED 
FREQUENCY  AND  DAMPING  RATIO  FOR  STRESS  ESTIMATION 
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APPENDIX  4A 
COMPUTER  PROGRAM 


4A.1  General  Notes 

The  fundamental  natural  frequency  and  r.m.s,  atreaaes  of  a  panel  with 
flexible  atiffenera  may  be  found  using  thle  computer  program.  The  locations  on  the 
plate  and  atiffener  at  which  r.m.a,  atreaaea  are  eatimated  are  given  in  Sketch  (ii). 
The  aasumptiona  given  in  paragraph  4,2  are  applicable  to  this  program. 

A  Hating  of  inatruotiona  for  two  sub-programa  is  given  in  FORTRAN  IV 
programming  language.  A  main  program  ia  required  to  read  in  data  and  print  out 
calculated  valuea  of  natural  frequency  and  r.m.a.  atresa.  A  listing  of  instructions 
ia  not  given  for  the  main  program  as  the  inatructions  required  are  dependent  on 
the  particular  computer  used.  Ths  requirements  for  the  main  program  are  given 
below. 


The  main  program  must  include  the  COMMON  statement  which  ia  written  in  both 
the  frequency  and  stress  sub-programs. 

4A.2  Frequency  Sub-Program 

This  sub-program  solves  for  the  fundamental  natural  frequency  of  the  plates. 
The  plate  edges  are  assumed  to  fixed  against  rotation. 

The  panel  frequency  data  to  be  input  for  each  panel  considered  are  values  for 
the  variables  listed  in  Table  4A.1. 


TABLE  4A.1 


Variable 

Variable  name 

Variable 

Variable  name 

A 

• 

AS 

I 

STI 

A 

A 

t 

T 

b 

B 

P 

RHO 

E 

EP 

p 

RHOS 

P 

• 

E 

ES 

8 

Any  coherent  set  of  units  in  which  time  is  expressed  in  seconds  may  be  used, 
the  frequency  being  obtained  in  Hs, 

On  returning  to  the  main  program  from  the  frequency  sub-program  the  fundamental 
natural  frequency  is  stored  in  variable  FREIJ . 

4A.3  Stress  Sub-Program 

This  sub-program  solves  for  the  r.m.s.  stress  in  the  plate  and  stiffener.  The 
geosietric  and  physical  property  data  required  are  the  data  for  the  frequency 
sub-program  with  the  additions]  variable  d  in  variable  name  DST.  The  panel  stress 
data  to  be  input  for  each  etrese  case  considered  are  values  for  the  variables 
listed  in  Tables  4A.3. 


TABLE  4A.3 


Variable 

Variable  name 

L  (f) 

SPL 

ps' 

£> 

DELTA 

In  this  sub-program  the  r.m.s.  fluctuating  pressure  equivalent  to  L  (f)  is 

ps 

2 

computed.  The  units  of  this  pressure  in  the  sub-program  are  N/m  .  If  British  units 

2 

are  used  this  pressure  must  be  calculated  in  units  of  lbf/in  . 


To  obtain  this 


pressure  in  Briti.h  unit.  th.  .ft— »«  against  l.b.l  106  .hould  b.  r.pl.c.d  by 

Y2*10.0*» (SPL/20 .0-8. 53749) 

On  returning  to  th.  main  program  from  the  atr.a.  aub-progra*  th.  r.m... 
•tresses  .r.  stored  in  ARRAY  STR  a.  shown  in  Tab.l  4A.4. 

TABLE  4A.it 


Streaa  location 

r- 

1 

2 

3 

4 

5 

Store  location 

str(i) 

STR ( 2 ) 

STR(3) 

STR(4 ) 

STR(5) 

o  a  o  |c«  a  n  o 


5: 


SUBROUTINE  FOR  PANEL  REFERENCE  FREQUENC  ' 

SUBROUTINE  FREQY 

COMMON  AS,A,B,EP,ES,STI,T,RHO,RHOS,DSTfSPL,DELTAfFREQ,STR(6),D,SK 

PNU-0.3 


CALCULATE  PLATE  REFERENCE  FREQUENCY 

ALPHA-3 ,0« { A/B ) **4+2 .0* ( A/B )  **  2 
D»£P*T**3/(l2.0*(l.O-PNU**2)) 

BETA-RHOS*  AS/ (RHO*  B*T)-ES*STI/ (D  *B ) 
GAMHA»AL?HA*(1.0+RH0S*AS/(RH0*B*T) )+3.0*BETA 
XI -SQHT( (GAMMA/ (2,0*  ALPHA ) ) **  2- 2 . 0*BETA/ALPHA ) 
X2— GAMMA/  ( 2 . 0* ALPHA  )  +X1 
X3— GAMMA/  (  2  .0* ALPHA  j  -XI 
X4=(X2*(ALPHAO.0)+2.0)/(3.0*X2+2.0) 
X5»(X3*(ALPHA+3.0)+2.0)/(3.0*X3+2.0) 
IF(X5)100,100,101 

101  IF(X4-X5)l00,100,102 
100  SKF-StlRT(X4) 

SK-X2 
GO  TO  103 

102  SKF=SQRT(X5) 

SK-X3 

103  FREQ-6 . 28 3 2*  S KF* ( S QRT ( D/ ( 3 . 0*  A **  4*  RH0*T ) ) ) 
RETURN 

END 


UBROUTINE  FOR  PLATE  AND  SIFIENEH  R.M.S.  STRESS 


SUBROUTINE  STRESS 

COMMON  AS , A  B  EP , ES , STI , T , RHO , RHOS , DST , S PL , DELTA , FREy , STR ( 6 ) , D , S  K 
DIMENSION  X(6) 

PNU=0 , 3 


CALCULATE  RMS  STRESSES 


Y=EP*T*A**  2*(l  ,0+SK)/(l9 .739^*  (l  .O-PNU  *'■  2 )* D  *(SK*4  2*  (  3 .0*  (  A/B )  > 
12,0*  (A/B)**  2  +  3 ,0)+4.0*SK+2,0*(l,  O+ES*  STI/  ( D*B)  )  )  ) 

X(l )=2.0*SK*(a/b)**2-PNU 
X(2)=1.0-2.0*SK*(A/B)*+2* PNU 
X(3)=1.0+2.0«SK 
X(4)=X(3)+2.0*PNU*SK* (  A/  H )  *  *  2 
X(5)=2.0*S»f  (a/ B  )  **  2  +  PNU*  ( 1 , 04  2 . 0*SK ) 

X(6)=ES*DST*(1 ,0-PNU**  2 )/ (EP*T ) 

IF(DELTA)104,I04,105 

104  DELTA=0 , 017 

105  Yl  =  (  0.78  54/DELTA)  *-*0.5 

106  Y2=10,0+*(SPL/20, 0-4. 69897) 

Y3-Y1* Y2*  FHEl}**  0 . 5*Y 

DO  107  1=1 ,6 

107  STH(l)=Y3*ABS(X(l)) 


RETURN 


END 


4  + 


